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ABSTRACT 


In  Part  II  of  this  paper,  calcvilatlons  are  made  for  the  bottom  loss 
of  a  plane  sound  wave  reflecting  from  a  layered  model  of  the  ocean 
bottom.  The  sea  floor  is  considered  to  consist  of  an  arbitrary  number 
of  absorbing  solid  layers,  each  of  constant  thickness,  lying  on  an 
absorbing  solid  half-space. 

Ccoqparisons  are  made  between  the  calculated  losses  and  the  experi¬ 
mental  values  given  in  Part  I.  In  one  of  the  experimental  areas  a  core  saa^Jle 
of  the  sediment  was  obtained  so  that  a  bottom  model  could  be  constructed 
to  accurately  represent  the  ^  situ  layering.  Attenuation  constants  for 
this  model  were  chosen  so  that  the  calculated  and  experimental  bottom 
loss  curves  are  in  reasonable  agreement.  These  constants  are  in  general 
agreement,  assuming  a  linear  relation  between  sound  attenuation  (in  db/ft) 
and  frequency,  with  high  frequency  resonant  chaotber  measurements  on 
similar  sediment  types.  In  the  other  two  areas  less  is  known  about  the 
natiire  of  the  sea  floor  sediments.  Bottom  models  constructed  for  these 
areas  provide  possible  explanations  of  the  experimental  loss  curves. 

A  short  discussion  of  bottom  loss  fluctuations  concludes  Part  II. 

A  study  of  the  possible  sound  paths  leads  to  the  conclusion  that  most  of 
the  fluctuations  observed  in  two  of  the  three  areas  are  probably  due  to 
Interference  between  the  four  possible  signal  paths  having  a  single 
bottom  reflection.  These  fluctuations  should  be  expected  in  any  test 
where  the  pulse , lengths  and  source  and  receiver  depths  are  such  that,  the 


A-\ 


bottom  reflected  signals  are  not  resolved.  At  the  third  area  the 
bottom  loss  fluctyations  are  large.  This  is  most probably  the  result 
of  combining  experimental  measurements  from  slightly  different 
locations  in  an  area  where  there  are  strong,  but  somewhat  irregular, 
reflecting  surfaces. 
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INTRODUCTION 

If  a  Bound  wave  is  projected  toward  the  ocean  bottom  there  will  be  a 
specularly  reflected  wave  of  diminished  amplitude.  The  ratio  i( in  db)  of 
the  soxjnd  intensity  of  this  reflected  wave  to  that  of  the  incident  wave 
is  called  the  bottom  loss. 

Part  II  of  this  paper  is  an  attempt  to  explain  the  experimental 
bottom  losses  reported  in  Part  I,  in  terms  of  available  information  about 
the  test  conditions  and  the  bottom  sediments  present  in  the  three  areas. 

For  this  purpose,  the  first  and  largest  section  of  Part  II  is  devoted  to 
the  development  of  equations  which  can  be  solved  for  the  bottom  loss  of 
a  plane  sinusoidal  wave  in  an  isotropic  liquid  half-space  and  incident 
upon  a  number  of  isotropic  absorbing  solid  layers,  each  of  consteuit 
thickness,  lying  on  an  isotropic  absorbing  solid  half-space.  T^iis  sample 
model  has  obvious  limitations  in  its  representation  of  any  real  bottom 
reflection  problem,  however,  in  most  cases  it  wou.:*!  serve  as  an  approximate 
solution  to  which  the  effects  of  scattering,  deviations  due  to  non-plane 
wave  front!,  non-sinusoidal  signals,  uneven  sediment  layering,  etc.  can 
be  addled. 

Following  the  development  of  equations  for  the  bottom  loss,  calcula¬ 
tions  are  made  for  the  areas  discussed  in  the  second  section  of  Part  I. 
Compaxison  of  calculated  and  experimental  values  are  presented  in 
graphical  form. 
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In  the  third  section  a  brief  discussion  of  the  fluctuation  of  the 
bottom  loss  values  is  given.  The  ocperimental  fluctuations  measured 
in  the  deep  water  areas  (II  and  III)  are  coinpared  with  the  expected 
four-signal  fluctuation,  derived  in  the  Appendix. 

EQUATIONS  FOR  THE  BOTTOM  LOSS 

The  ocean  bottom  has  been  traditionally  treated  as  a  liquid  in 

1  2 

problems  invol.ing  bottom  reflection  '  .  When  the  absorption  of 
sound  in  the  bottom  needs  to  be  considered,  and  when  this  absorption 

is  small,  a  convenient  model  to  use  is  that  of  a  liquid  with  a  complex 

3  4 

velocity  of  sound  .  In  a  recent  paper,  T.  G.  Bell  reported  on  the 

comparison  of  experimental  bottom  losses  in  several  areas  with  theoretical 

values  calculated  using  a  two-layer  "absorbing"  liquid  model.  The 

agreement  between  the  calculated  and  observed  values  of  bottom  loss  is 

excellent.  However,  the  neglect  in  these  studies  of  the  shear  waves 

that  would  be  generated  in  the  viscous  liquid  layers  might  be  questioned, 

since  in  several  of  the  models  the  attenuation  of  the  compressional  wave 

was  of  the  order  of  0.3  -  0.6  db/ft.  A  more  realistic  model  is  that 

of  a  visco-elastic  solid^  in  which  the  effects  of  bottom  rigidity  and 

absorption  are  in  c luded. 

The  problem  of  calculating  the  reflection  coefficient  of  a  sound 
wave  incident  on  a  visco-elastic  layer  is  essentially  the  same  as  for 
a  wave  incident  on  an  elastic  layer^,  except  that  the  Lam6  constants 
are  complex  numbers. 
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It  Is  assumed  here  that  an  area  of  ocean  bottom  ^diich  is  under 
consideration  can  be  represented  as  a  set  of  isotropic,  constant 
thickness,  visco-elastic  solid  l£^yers.  Hie  bottom  layer  has  infinite 
thickness,  i.e,  it  is  a  half-space.  A  plane  sound  wave  in  the  water 
layer  and  incident  on  this  model  bottom  will  generate  in  each  layer, 
with  the  exception  of  the  lowest  layer,  four  sets  of  plane  waves^  two 
sets  of  up-and-down  going  coiapressional  Waves  and  two  similar  sets  of 
shear  waves.  In  the  bottom  layer  only  the  down-going  waves  are  present. 

As  in  Figure  1,  denote  the  n^^  layer  by  the  same  subscript.  The 
sets  of  waves  in  the  n  layer  can  be  defined,  in  the  usual  manner,  in 
terms  of  four  scalar  velocity  potential  functions'^: 

P-  - 

n  '  • 

- 1  , 

'Pn^Ap  exp.  ;  ( kx 4- Op kz -  a;t ) , 

^n=An  exp.  ((kx-Opkz-'cot), 
v//^  .=  Bp  exp.  i  (kx +bn  kz  -  ojt )  1 
\;/p  =  Bp  exp.  i  (kx  -  bpkz  -coty  .' 

It  is  assumed  that  only  real  values  of  k  need  be  considered  to 

obtain  a  good  approximation  to  the  incident  wave  in  the*  water.  In  this 

case  k  is  determined  by  the  angular  frequency  ou,  the  angle  of  incidence 

of  the  incoming  sound  wave  0  and  the  velocity  of  sovind  in  the  water  c 

1»  *'  w  . 

by  the  relation  k  =  o)  sin  0,/c  .  The  values  of  a  and  b  are  determined 

1  w  n  n  . 


by  the  physical  properties  of  the  layer  and  the  condition  that  the 
potentials  must  satisfy  their  respective  wave  equations.  It  follows 
that: 

“n  =[  +  -  l> 

]-l  ■ 

Here,  is  the  density  of  the  n^^  layer  and  c,  the  horizontal  canponent 

of  the  phase  velocity,  Is  equal  to  oy^k.  The  vlsco-elastlc  properties 

of  the  layer  are  determined  by  the  complex  Lame  constants  \  and 

Cotislstertt  with  the  assumed  time  dependence  of  the  potential  functions, 

X  ceui  be  written  as  -  1  X"  and  u,  as  u'  -  1  u"  where 

n  n  n  V.  n  '^n  n 

2 

X^,  X|^,  and  are  positive  real  numbers.  Ihe  roots  of  a^ 

2 

and  b^  should  be  chosen  to  lie  in  the  first  quadrant  of  the  complex 

pl6me.  In  order  to  reduce  the  niomber  of  constants  required  to  define 

the  model  it  will  be  assianed  here  that  the  volume ,  viscosity  is  equal 

to  zero,  that  is,  X"  =  -2m,”  /  3. 

'  '  n  n  ' 

At  the  interface  between  the  two  visco-elastic  layers  the 
components  of  particle  velocity  in  the  vertical  (z)  and  horizontal  (x) 
directions,  respectively  denoted  by  w  and  u,  and  the  elements  p 

zz 

and  p  of  the  stress  tensor  must  be  continuous.  In  the  case  of 
xz 

plane  sinusoidal  waves,  assumed  here,  the  relations  between  the 

canponent  velocities,  the  stress  tensor  elements  and  the  velocity 

« 

potentials  can  be  written  as: 
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Referring  to  Fig.  1  let  the  n  interface  be  a  distance 
the  water-sediment  interface.  It  follows  that  the  interface 
at  the  n^^  interface  (at  depth  z*!h)  can  be  written  as; 


Op  [  exp.  i  (apkh)-  a’^  exp.  i  (-a^kh)] 

exp.  i  (  bpkh  )  +  B|I|  exp.  i  (-b^kh  ) 

=  similar  expression  for  layer  (n  +  l). 

A'r,  exp.  i  (  a^kh  )  +  A^  exp.  i  (-  a^,  kh  ) 

■^n  [  Bn  exp.  i  (  bnkh)  -  B'n  exp.  i  (-  kh)  ] 
^similar  expression  for^layer  {n-fl  1. 

^n  ®*P  '  A'n  exp.  i(-a 

+  2  b„  [Bn  exp.  i  (bnkh)  -  Bp  exp.  i  {-bpkh 
=  similar  •  expression  for  layer  (n+l). 

/Ip  |-2ap[  Ap  exp.  I  (  apkh  )-  Ap  exp.  i  (-a, 
+  (bp  -  I  )  [  Bp  exp.  i  (bpkh)  Bp  exp.  i  (-bp 
s  simiar  expression  for  layer  (n+l) 


h  below 
conditions 


nkh)] 
>1} 

,kh)] 

kh)]}' 


•Hie  equations  valid  at  the  other  Interfaces  between  layers  can 
be  obtained  by  appropriate  change  of  subscripts  amd  the  depth  of  the 


Interface.  In  the  equations  associated  with  the  deepest  interface  the 
coefficients  and  are  equal  to  zero.  At  the  water-sediment  inter¬ 
face  (interface  p  at  z  c  o),  there  are  o:ily  three  Interface  equations 

V.  '  *  "p‘  ‘ 

'  ‘/^p  C  'CAp+Ali,)  . 

+2bp(B;,-B;i)]. 

2ap(Ap-Aj)»  (b|-l)  (Bp+  Bp), 


Here,  is  the  density  of  the  water  and  a  ,  can  be  Identified  as 

P+J-  p+1 

cot  6.  and  k  as  o)  sin 

1  .  1  ^  • 

In  general  the  Interface  conditions  will  result  in  Un  +  3  linear 
algebraic  equations  (N  =  number  of  layers)  relating  the  4n  +  4  complex 
amplitude  factors.  By  dividing  each  equation  with  A^^  the  number  of 
unknown  quantities  is  reduced  to  the  number  of  equations  and  a  solution 
can  be  obtained  for  the  reflection  coefficient  A"  ,/A'  .  The  bottom 

P+1'  p+1 

loss  in  db  is  minus  twenty  times  the  log  (base  lO)  of  the  modulus  of 
the  reflection  coefficient.  One  direct  method  of  solution  for  the 
reflection  coefficient  is  to  expand  the  set  of  con^jlex  equations  into 
twice  that  number  of  real  algebraic  equations  and  to  solve  for  the 


ref loc cion  coefficient  by  an  iterative  inversion  of  the  matrix  of  real 
coefficients.  However,  when  the  number  of  layers  becomes  large  the 
computation  time  required  to  set  up  and  invert  the  coefficient  matrix 
becomes  excessive  (the  time  required  to  inArert  a  matrix  in  most 
Iterative  methods  increases  as  the  third  power  of  the  order  of  the  matrix). 
An  efficient  method  of  calculating  the  reflection  coefficient  for  a  model 
with  a  large  number  of  layers  is  developed  in  the  following  paragraphs. 

From  physical  arguments,  the  stresses  and  particle  velocities  at 
one  interface  of  a  layer  are  related  to  the  stresses  and  velocities  at 
the  second  Interface.  For  small  displacements  these  relations  are 
linear  and  can'be  expressed  as  a  matrix  equation.  This  method  is  due 

Q 

to  W.  T.  Thomson  who  calculated  the  reflecti6n  coefficient  of  a  plane 

wave  in  a  liquid  incident  on  a  set  of  elastic  solid  layers  bounded 

below  by  a  liquid  half-space.  Here  the  same  procedure  is  followed,  but 

each  layer  and  the  lower  half-space  are  assumed  toba  visco-elastic 

solids.  PC  30,  a  minor  correction  to  Thomson's  equations  concerning 

the  interface  conditions  has  been  made. 

m  th 

Let  u^  be  the  x  component  of  particle  velocity  in  the  n  layer 

th 

and  at  the  m  interface .  Here  and  in  later  parts  of  this  paper  sub¬ 
scripts  will  refer  to  a  particular  layer  and  superscripts  to  a  specific 
interface.  Let  the  origin  of  the  coordinate  system  be  shifted  to  the 

n-1  interface  so  that  the  n^^  interface  is  located  at  z'*«  -d  ,  where  d 

n  n 

is  the  thickness  of  the  n^^  layer.  From  previous  equations  the  expression 
of  u^,  with  the  time  factor  emitted,  can  be  written  as: 


c  c 


=  lk([Anexp.i(-ankndn)  +  A';,exp.l  (ank^d^)] 

-bn  [Bn  ®xp.i  (-bnkndn)-BS  exp.i*  (bnkndn)]J  exp.l(kx). 

Similar  equations  can  be  written  for  the  z  con5)onent  of  particle 

velocity  in  the  n^^  layer  emd  at  the  n^^  interface,  and  for  the  stress 

tensor  elements  and  (p  .  Let  V™  be  a  column  matrix  with 

'  zz'n  xz'n  n 

elements  u“,  w“,  ^^xz^n^  vhere  m  may  be  n  or  n-1. 

Let  W  be  a  column  matrix  with  elements  A'  +  A",  A'  -  A", 
n  n  n^  n  n^ 

B'  -  B"  and  B'  +  B"  .  llien  V™  and  W  can  be  connected  by  a  matrix 
n  '  n  n  n  n  n 

.equation;  let  .  W  .  The  4x4  matrix  is  given  in  Table  1. 

n  Ti  n  n 

Because  the  interface  n-1  is  located  at  z  =  0  the  elements  of  can 

n 

be  obtained  frcm  the  elements  of  by  setting  d^  e  0.  Note  that  half 

of  the  elements  of  are  zero  because  of  the  plus  and  minus  arrangement 

of  the  A’s  and  B’s  in  W  .  The  column  matrices  and  are 

n  n  n 

related  by  the  equation  h  D  .  .v”~^  .  The 

n  n  n  n  ■  n  n 

elements  of  D  are  listed  in  Table  8. 
n 

Prom  physical  arguments,  the  relations  between  and  must 

n  n 

be  Independent  of  the  location  of  the  origin  of  the  coordinate  system. 

Therefore,  the  quantity  d  ,  which  occurs  in  the  elements  of  D  ,  has 

n  n 
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significance  only  as  the  thickness  of  layer  n.  Thus,  matrices  of 

the  same  form  can  he  written  for  the  other  layers.  The  boundary 

condition  at  Interface  n  can  be  expressed  as  It  follows 

n+l  n 


that  o  D  . 

P  P  P 

or  finally,  =  D 
P  p 


For  convenience,  the 


origin  of  the  coordinate  system  is  now  shifted  to  the  Interface  1.  Then, 
“i-  The  matrix  has  the  same  form  as  with  d^  «•  0 , 


because  sxirface  1  is  now  located  at  z  =  0.  For  convenience,  define 


P 


D  .  D  , 
P  P-1 


.  .  .  Dg  • 


s  F 


Let  the  total  thickness  of  the  layers  2  through  p  be  H.  The  water- 
solid  interface  is,  therefore,  located  at  z  =  -H.  Noting  thatihe  p-1 
half-speice  is  a  liquid,  the  following  interface  conditions  hold  at  the 
surface  p:  ■  . 


+  ®XP. •("Op^jkH) “Ap^l  exp, I (Op^jkH)J  =  Wp  , 

"‘'"p  +  l  "[Ap  +  i  e’<P.M-Op^,kH)+A‘;>^.,exp.l(Gp^,kH)^^^^^ 
^  '  ^Pxz^p  • 

( 

Let  R  =  \  =  V^'p+1^  ^s  "  ®i/'^'p+l*  preceding 

equations  can  then  be  expanded  Into  the  following  form  in  terms  of  the 


elements  of  the  matrix  F. 


ip+,kH)] 


iOp^jk  e)  i(-ap^|kH)[l  -  R  exp.i' 
■  (  ^21  ‘*’•^22^  ^^23^  ^24^  ^8  » 


"'^P+l  ^  R  (20p^jkH)] 


^  P^3J  '*’  ^Z2  ^  **"  ^  **'  ^  "^8  * 


33  ’  '  34'  *8 


0-  (F^,+  F;2)Tc  +  (F43  +  F44)Ts. 

This  Bet  of  equations  is  solved  to  find  R  exp.i(2a^^^kH),  which  is 

tl  I 

the  desired  ratio  of  reflected-to-incident  sound  pressure 

at  the  interface  between  the  water  and  the  top  sediment  layer  (at 
z  B  -H) .  The  bottom  loss  is  defined  as  minus  twenty  times  the  log 
(base  lO)  of  the  modulus  of  the  reflection  coefficient  R  exp  i(2a^^kH). 
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CCMPARISON  OF  EXPERIMENTAL  AND  CALCULATED 
BOTTOM  LOSSES 

Area  I,  discussed  in  Part  I  of  this  paper,  lies  in  moderately 
deep  water  off  the  coast  of  Southern  California.  The  physical  properties 
of  the  bottom  in  Area  I  Jiave  been  studied  by  a  group  from  the 
U.  S.  Navy  Electronics  Laboratory  under  E..L.  Hamilton.  Low  frequency 

Q 

sub-botlom  profile  recordings  by  D.  G.  Moore  indicate  a  strongly 
reflecting  layer  approximately  20  feet  below  the  water-bottom  interface. 
As  noted  in  Part  I  of  this  paper,  the  ccoipositlon  of  the  upper  sediment 
layers  varies  throughout  the  operating  area.  The  bottom  model  assumed 
in  this  area  is  based  on  several  core  samples,  obtained  by  Hamilton. 

The  coii5)ressional  velocities  in  the  layers  fovind  in  the  core  sample^  were 
determined  by  the  resonance  method  of  Shizmway^*^''  ^  and  are  listed  in 
column  9,  Ifeible  3#  Model  A. 

A  mathematical  model  of  the  bottom  was  then  constructed  (Table  3, 
Model  a)  for  Area  I.  The  experimental  values  of  the  velocity  and 
density  of  a  layer  essentially  determine  the  value  of  (X'  +  2pi')  for 
the  layer^.  A  value  was  then  chosen  for  the  ratio  m^'A'  according 
to  the  type  of  material  found  in  the  core  sampl^.  In  the  models 
described  in  this  paper,  the  choice  of  this  ratio  is  not  critical. 
Variations  in  the  value  of  (p,'A')  of  the  order  of  0.1  has  small  effect 
on  the  bottom  J.OSS  curves.  It  may  be  noted  that  choosing  a  value  for 


pi'/X'  iB  equivalent  to  specifying  Poisson's  ratio  for  the  sediment 

layer.  The  attenuation  constants  were  chosen  to  provide  reasonable 

agreement  between  the  experimental  and  calculated  values  of 

bottom  loss.  The  values  are  in  general  agreement,  assuming  a  linear 

relation  between  sound  absorption  (in  db/ft)  and  frequency,  with  the 

high  frequency  resonant*  chamber  measurements  of  Shumway^^  on  similar 

types  of  sediment.  In  contrast,  however,  these  attenuation  constants 

resvilt  in  a  much  higher  attenuation  of  the  compressioneil  wave  than 

12 

that  measured  by  Wood  and  Weston  in  mud.  The  bottom  model  constants 
(X',  Vi',  Vi",  and  p)  are  listed  in  Table  3,  Model  A.  A  comparison  of 
the  calculated  and  experimental  bottom  losses,  as  a  function  of 
grazing  angle,  is  shown  in  Fig.  2.  Calculations  for  Fig.  2,  and  the 
other  figures  in  this  paper,  were  made  using  both  the  matrix  inversion 

and  the  jnatrix  multiplication  type  solutions  for  the  bottom  loss,  discussed 
earlier.  For  all  models  and  at  all  grazing  angles  the  difference  between 
the  calculated  values  of  bottom  loss  by  the  two  methods  of  solution  was 
less  than  0.01  db.  The  curved  line  represents  the  calculated  bottom 
losses  and  the  straight  lines  connect  points  which  are  average'  values  of 
the  experimental  losses  combined  in  twh-degree  intervals.  The  caps  on 
the  vertical  lines  represent  values  of  bottom  loss  equEil  to  the 
mean  experimental  value  plus  or  minus  the  standard  deviation. 

In  Fig.  3  a  compariscn  is  made  between  experimental  values, 
averaged  over  five-degree  intervals,  and  calculated  bottom  losses  for  a 
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two-layer  model,  the  upper  layer  six  feet  thick*.  This  model  is  listed 
as  Model  B  in  Table  3.  A  study  of  Figs.  2  and  3  led  to  the  following 
conclusions.  The  broad  peaks  shown  in  the  experimental  losses 
averaged  over  five-degree  intervals  are  the  result  of  interference  of 
sound  reflected  from  the  surface  of  the  bottom  and  from  the  interface 
at  approximately  six  fdet.  The  rapid  fluctuations  that  appear  in  the 
experimented,  curves  averaged  over  two-degree  intervals  are  probably  due 
to  the  interference  as  a  result  of  sound  reflecting  from  the  20-ft  deep 
interface.  The  large  standard  deviations  of  the  individual  data  points 
are  discussed  in  the  following  section. 

In  contrast  to  Area  I,  core  samples  were  not  taken  in  Areas  II  and 

III.  Consequently,  the  mathematical  models  of  the  sediment^  in  these 

S 

areas  are  overly  simple.  The  layer  constants  and  jx*  were  chosen  so 
that  the  calculated  and  experimental  values  of  bottom  loss  would  have 
the  some  general  angular  dependence  and  so  the  calculated  compressional 
velocity  in  the  layer  would  have  a  reasonable  value.  The  attenuation 
constants  were  then  adjusted  for  agreement  between  the  experimental  and 
calculated  bottom  loss  curves  at  small  grazing  angles. 

Area  II  is  deep  water  area  off  the  California  coast.  The  bottom  is 
believed  to  be  a  silty  clay.  Because  a  detailed  bottcmi  study  was  not 

t 

made  in  this  area,  a  single  layer  model  was  chosen  for  the  bottom  loss 

*Before  the  present  calculations  were  begun  and  before  core  samples  in 
this  area  were  available,  Mr.  Garland  Bernard  of  the  Defense  Research 
Laboratory,  Austin,  Texas  suggested  that,  on  the  basis  of  bottom  loss 
calculations  for  a  two-layer  absorbing  liquid  bottom,  the  experimental 
bottom  losses  were  consistent  with  a  discontinuity  in  the  bottom  at  a 
depth  of  six  feet.  Subsequent  core  samples  proved  this  to  be  true. 
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calcxilatlons.  Dm  experimental  tests  in  Area  II  were  made  at 

frequencies  of  O.J,  1.3>  &nd  3*0  kc.  Because  the  set  of  measurements 

at  1.3  kc  covered  a  larger  range  of  grazing  angles  the  model  constants 

were  adjusted  for  this  frequency.  Hieoretical  curves  at  0.7  and  3.0  kc 

were  then  calculated  assuming  two  different  frequency  variations  for 

the  attenuation  constant  (i",  as  explained  later.  The  experimental  data 

Indicate  rather  large  overall  hottcm  losses  with  a  maximum  bottom  loss 

for  the  1.3  kc  tests  at  a  grazing  angle  of  approximately  l4  degrees. 

The  maximum  loss  at  small  grazing  angles  can  be  explained  by  either  a 

low  velocity  bottom  (Model  C  and  Fig.  k)  or  a  highly  attenuating  bottom 

(Model  P  and  Fig.  5).  In  the  low  velocity  cases  the  high  bottom  loss 

Is  due  to  refraction  of  the  sound  into  the  bottom,  mostly  In  the  form 

of  a  compresslonal  wave,  while  in  the  high  attenuation  models  the  hi^ 

loss  Is  largely  due  to  conversion  of  sound  energy  into  a  shear  wave. 

That  widely  differing  models  can  be  used  to  explain  the  bottcm  loss 

curves  proves  the  need  for  measurements  of  the  physical  prc^rties  of 

the  bottom,  partlculaurly  the  ^  situ  ratio  of  the  soxind  velocity  in 

water  to  the  compresslonal  velocity  In  the  upper  few  feet  of  sediment. 

This  ^  situ  ratio  hsis  been  determined  in  a  few  cases  by  Hamlltcsi 

13 

using  the  bathyscaph  TEIIESTE  . 

For  both  the  low  velocity  models  and  the  high  velocity  models  the 
attenuation  constant  n"  was  chosen  for  a  fit  of  the  experimental 
and  calcxilated  bottom  loss  curves  at  1.3  kc.  The  theoretical  loss 
curves  at  0.7  and  3.0  kc  are  determined  by  the  frequency  dependence  of  4" 


It  seems  reasonable  that  over  this  frequency  range  can  be  assumed 
to  vary  as  n"  ■  equal  to  a  constant  and  OSn^l.  For  example, 

if  n  is  0  the  velocities  of  the  compress ioneil  and  shear  waves  are 
independent  of  frequency  and  the  attenuation  of  these  waves  will  be 
linear  with  frequency.  If  n  is  1,  the  sediment  layers  are  represented 
as  a  Voig  t  solid  type*  material  and  the  velocity  and  attenuation  rela¬ 
tions  are  more  ccnqpllcated^ , 

In  this  paper,  it  was  assumed  alternatively  that  and 

u"  ■  (n  ■  1  or  O).  Models  C,  D,  and  E  (Table  3  ancl  Fig.  4) 
represent  the  low  velocity  case,  with  ji"  =  at  frequencies  of  1.3/ 

0.7/  and  3*0  kc  respectively.  Note  that  Model  E  is  not  a  low  velocity 
case  since,  as'  a  result  of  the  increase  in  the  attenuation  constant, 
the  conpresslonal  velocity  in  the  sediment  is  higher  than  the  water 
velocity.  For  the  assun^tlon  u"  ■  the  losses  are  independent  of 

frequency  andNthe  bottom  loss  curves  are  given 'by  Model  C,  Fig.  4. 

Models  F,  0,  and  H  (Table  3  and  Fig.  5')  have  the  same  correspondence  for 
the  high  velocity  models. 

Itoe  experimental  losses,  as  shown  in  Fig.  4  or  5/  do  not  increeise 
at  the  higher  frequencies,  however,  it  cannot  be  concluded  that  the 
assuoption  n  =  1  is  a  better  choice.  At  least  for  many  oceaui  bottoms, 

t 

the  single  layer  model  is  a  poor  representation  of  the  sediment  layers, 
particularly  for  low  frequency  bottom  loss  calculations.  A  multi-  layer 
model  will  show  a  general  increase  of  bottom  loss  with  frequency  for 


either  assumption.  It  might  be  noted  that  the  assumption  n  >  1  leads 

to  Boae  unlikely  Increases  of  the  c empress tonal  velocity  with  frequency 

(Table  3*  Column  9)«  Also,  there  are  reports  on  the  direct  measurement 

of  attenuation  constants  which  indicate  that  the  attenuation  terms  are 
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Independent  of  fijequency  * 

Area  III  is  a  deep  water  area  of  the  Bering  Sea.  The  reflection 
measurements  immediately  sxiggest  that  the  low  losses  at  angles  near 
normal  Incidence  are  due  to  a  heu*d  reflecting  layer,  such  as  basalt. 

The  larger  losses  at  small  grazing  angles  are  most  likely  due  to  sm  ab¬ 
sorbing  layer,  such  as  silt,  lying  on  top  of  the  hard  layer.  A  two- 
layer  model  was  used  and  the  attenuation  of  the  upper  layer  was  adjusted 
to  fit  the  peak  of  the  experimental  curve  at  the  1.3  kc  frequency  (Fig.  6). 
In  Pig.  7  the  1.3  kc  curve  is  repeated  along  with  experimental  and 
calculated  curves  for  the  bottom  loss  at  0.7  kc.  The  0.7  kc  curves  were 
calculated  assvaning  ji"  to  be  constant  and  then  to  be  a  linear  function 
of  frequency.  The  agreement  with  experimental  values  is  only  fair  for 
either  evurve.  Experimental  data  at  3-0  kc  were  not  obtained  in  this  area. 

SOME  COMMENTS  ON  THE  FLUCTUATIONS  OF  THE  BOTTOM  LOSSES 
Fluctuations  in  the  Intensity  of  bottom  reflected  sound  are  usually 
attributed  to  interference  of  signals  traveling  along  different  paths,  or 
scattered  from  reflecting  objects  lying  on  the  bottom,'  or  scattered  from 
reflecting  facets  of  the  surface  of  the  bottom.  In  the  areas  discussed 


n 


in  this  paper,  the  water-sediment  interfaces  are  relatively  smooth  so 
that  it  was  expected  that  fluctuations  of  the  bottom  loss  would  be  small. 

Standard  deviations,  averaged  over  all  grazing  angles,  of  the  bottom 
loss  measurements  in  Area  II  were  3.8,  3.2  and  3.3  db  at  frequencies  of 
0.7,  1.3>  arid  3*0  kc  respectively.  In  Area  HI  the  standard  deviations 
of  the  bottom  loss  were  2.4  db  at  I.3  kc  and  1.6  db  at  0.7  kc.  In  the 
deep  water  areas  (ll  and  III)  shallow  source  and  receiver  depths  and 
long  pxilse  length  (0.5  sec)  signals  were  used.  The  fluctuations  result, 
most  probably  from  scattering  and  from  mutual  interference  of  the  four 
possible  signal  paths  having  a  single  bottom  reflection,  which  are 
described  in  the  Appendix.  In  moderate  seas  the  surface  motion  of  the 
water  will  cause  the  difference  in  path  lengths  of  the  four  signals  to 
vary  by  more  than  a  wavelength  at  these  frequencies.  Thus,  between  suc¬ 
cessive  pings  the  phase  differences  of  the  foxcc  signals  will  change  in  a 
random  manner.  -  . 

Results  from  the  Appendix  show  th^t  a  standard  deviation  of  about 
2.3  db  is  expected  because  of  four-path  interference  of  constant  amplitude 
signers  in  Areas  II  and  III.  Since  each  of  the  four  signals  will  have 

y 

some  amplitude  fluctuation  the  actual  four-path  inierference  should  be 
somewhat  larger.  •  Thus  it  seems  that  in  Area  II  a  large  part  of  the 

■  t 

observed  fluctuations  is  due  to  four-signal  interference  and  in  Area  III 
almost  all  of  the  variation  in  bottom  loss  was  due  to  this  effect.  The 
low  value  of  fluctuation  (I.6  db)  at  0.7  kc  in  Area  III  can  be  explained  by 
noting  that  these  tests  were  made  in  unusually  calm  seas  and  that  at  the 


0.7  kc  frequency  the  phase  differences  hetven  the  interfering  signals 
were  not  randan. 

Ihe  experimental  losses  in  Area  1  were  measured  in  l4  different 
tests  made  over  a  period  of  several  months.  Thus  the  losses  reported 
for  a  certain  grazing  angle  were  obtained  at  several  locations  in  the 
same  general  area.  Because  of  this,  and  apparently  because  of  the 
strong  interference  of  sound  reflecting  from  different  surfaces  of  the 
sediment  layers,  the  fluctuations  of  the  bottan  loss  were  large  in  this 
area.  The  standard  deviation  was  4.4  db  when  the  losses  were  grouped 
in  two-degree  Intervals  and  4.9  db  when  taken  in  five-degree  Intervals. 
Smaller  samples  of  data  also  showed  considerable  fluctuation. 

CONCLUSION 

Equations  have  been  derived,  in  two  different  forms,  for  the  cal¬ 
culation  of  the  bottom  loss  of  a  plane  sound  wave  In  water  and  incident 

on  an  arbitrary  number  of  visco-elastic  solid  layers.  Thus  the  model 

of 

Includes  the  effects  sediment  layering,  the  absorption  of  sound  waves 
in  the  bottan  sediments,  and  the  generation  of  shear  waves  at  the 
Interfaces.  In  most  of  the  models  considered,  the  bottan  losses 
were  high  at  small  grazing  angles  because  of  a  Isurge  part  of  the 
so\md  energy  being  converted  into  a  shear  wave. 


The  botton  loss  as  a  function  of  grazing  angle  Is  calculated  for 
two  locations  off  the  California  coast  and  at  a  deep  water  location  in 
the  Bering  Sea.  In  only  one  of  these  locations  was  more  than  half  of 
the  bottom  constants,  required  to  define  the  model  bottom,  available. 
Ihus,  no  claim  can  be  made  about  how  well  this  model  will  represent  an 
ocean  bottom,  even  in  the  special  case  where  the  interfaces  are  smooth 
and  the  layering  uniform.  At  one  location  )C  core  seimple^'ww  available 
so  that  the  most  critical  assumptions  were  the  values  of  the  layer 
attenuation  constants.  The  chosen  values  of  these  constants  were  in 
general  agreement,  assuming  a  linear  relation  of  attenuation  to  frequency 
to  those  determined  at  high  frequencies  for  similar  sediments  by  the 
resonance  method. 

Fluctuations  in  the  bottom  loss  at  the  deep  water  locations  are  moat 
readily  explained  as  resulting  from  interference  of  the  four  signals 
with  one  bottom  reflection.  The  large  fluctuation  in  the  bottom  loss 
in  the  moderate  depth  water  location  is  most  probably  the  resxilt  of 
combining  experimental  measurements  from  slightly  displaced  stations  in 
an  area  where  thei*e  are  strong  but  somewhat  irregxilar  scattering  siirfaces 
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Appendix 


In  the  deep  water  areas  (ll  and  III)  discussed  in  this  report, 
the  source  and  receiver  were  both  located  at  shallow  depths  and  0.5  sec 
pxilses  were  used.  Under  these  conditions  the  four  signals  that  have  one 
bottcm  reflection  will  overlap.  Let  the  signal  that  is  reflected  from 
only  the  bottom  be  labeled  as  B,  the  signal  reflected  from  the  surface 
and  then  the  bottom  be  called  S-B,  etc.  Because  the  source  and 
receiver  were  located  at  nearly  eq.ual  depths  the  S-^B  and  B-S  paths  have 
almost  equal  travel  times.  Then  the  signal  recorded,  diuring  the  arrival 
of  the  four  single  bottan  reflected  paths  will  at  increasing  times  be 
made  up  of  the  following  five  combinations; 

1.  B 

2.  B,  S-B,  B-S 

3.  B,  S-B,  B-S,  S-B-S 

4.  S-B,  B-S,  S-B-S 

5.  S-B-S. 

In  the  analysis  of  the  experimental  records  the  peak  value  of  the 
recorded  signal  was  measured  which  may  occur  in  any  of  five  possible  com¬ 
binations  of  signals.  Let  f^(p)dp  be  the  probability  of  finding  a 
value  of  peak  pressure  between  the  values  of  p  to  •  p+dp  in  the  i^^ 
combination  of -signals.  Then  the  probability  of  finding  the  largest 
value  of  peak  pressure  from  any  of  the  five  combinations  to  be  in  the 

t. 

p  to  p+dp  is 


i^  ^  Jol/i  J  ^ 


y. 


It  is  assumed  now  that  each  of  the  four  interfering  Signals  has 


constant  aanplitude  and  that  the  phase  differences  between  the  interfering 


signals  ax-e  random.  The  distribution  functions  f^  and  f^  are  Dirac  . 
delta-functions.  The  distribution  functions  for  the  three  Interfering 


signals  (fg  and  f^)  can  be  written  as  complete  elliptic  Integrals  of 

* 

the  first  type . 


Let  the  B,  S-B,  B-S,  and  S-B-S  signals  have  the  respective  peak 

2 

amplitudes  .of  1,  a,  a,  and  ,  where  a^i.  Then  fg  is  the  distribution 

I 

function  of  the  ptak  amplitudes  1,  u,  and  a  combining  with  random  phase. 
Pig.  8  is  a  graph  of  fg  for  values  of  the  parameter  a  equal  to  1.0, 
0.8,  and  0.6. 


The  function  f^,  which  represents  the  combination  of  all  four  signals, 

17 

was  calculated,  following  the  method  of  Pearson  ,  by  numerical  integra¬ 
tion  of  a  distribution  function  for  three  signals.  Pig.  9  is  a  graph  of 
f^  for  values  of  a  equal  to  1.0,  0.8,  and  0.6.  *  The  functions  f^, 

and  f^^  can  also  be  calculated  in  terms  of  an  infinite  integral  or  as  an 
infinite  series  .  These  methods  were  not  used  here  because  of  their  slow 


convergence  when  the  number  of  combining  signals  is  small. 

The  desired  function  P^  is  easily  calculated  when  f^,  f^,  and  fj^ 

are  known.  Pig.  10  is  a  graph  of  P^  for  values  of  fl  equal  to  1.0,  0.8, 
and  0.6.  From  F^  the  distribution  function  of  the  bottom  loss  can  be 
calculated.  Although  the  value  of  the  bottom  loss  (in  db)  depends  upon  the 
reference  intensity  the  standard  deviation  of  the  bottom  loss  does  not.  The 


expected  standard  deviation  of  the  bottom  loss,  for  the  above  assvimptlons. 
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